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We introduce a quantum channel describing photon losses depending on the degree of polarization 
mixedness. This can be regarded as a model of quantum channel with correlated errors between 
discrete and continuous degrees of freedom. We consider classical information over a continuous 
alphabet encoded on weak coherent states as well as classical information over a discrete alphabet 
encoded on single photons using dual rail representation. In both cases we study the one-shot 
capacity of the channel and its behaviour in terms of correlation between losses and polarization 
mixedness. 

o\ " 

O 1 PACS numbers: 03.67.Hk, 42.50.Ex 

O 

<N 

; I. INTRODUCTION 

o . 

In recent years lots of attention has been devoted to capacity of quantum channels and related properties |l[ . Along 
this direction introducing models with correlated noise has provided more realistic framework for studying quantum 
channels (2j. Correlated noise was first introduced in [3[ for depolarizing channel, then in subsequent works [3, HJ 
the advantage of using entangled states was addressed for different channels and some of the effects has also been 
experimentally demonstrated [6(. The notion of quantum channels with correlated noise has also been extended to 
continuous- variable, e.g. lossy bosonic channels 7] and channels with additive noise [|[. 

In all the earlier works, correlation is defined between different uses of the channel, but within the same degree of 
freedom. In practice correlated errors may even arise in a single channel use between different degrees of freedom. 
Optical fibers, which are among the most successful candidates for quantum communication, are examples of this 
kind. In any use of an optical fiber, different types of noise may happen like loosing photons to the environment 
t-H ■ and at the same time having polarization effects. Loss effects in optical fibers can be modeled by a beam splitter 
i> interaction with a vacuum (or thermal) environment. The classical capacity of the memoryless lossy bosonic channel 
was exactly evaluated in Q, while for the same channel with memory the evaluation turns out to be more involved 
and to be model dependent 0- 

Although the capacity of lossy bosonic channels has been widely investigated, the study of polarization effects 
in the rate of information transmission is more restricted to classical optical communication [Toj . In this paper we 
introduce a quantum channel describing photon losses depending on the degree of polarization mixedness. This can 
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be regarded as a model of quantum channel with correlated errors between different degrees of freedom. We consider 
classical information over a continuous alphabet encoded on weak coherent states as well as classical information over 
a discrete alphabet encoded on single photons using dual rail representation. In both cases we study the one-shot 
capacity of the channel and its behaviour in terms of correlation between losses and polarization mixedness. 

The layout of the paper is the following. In section [TT| we introduce the model. We then study in section IIIII the 
behaviour of Holevo information versus the correlation parameter for the case of information encoded in weak coherent 
states. In contrast, in section HVl we use the dual-rail encoding for qubits. In this case our model realizes a quantum 
erasure channel [ll| where the degree of correlation controls the probability of erasing information. Paper ends with 
conclusion in section IVl 



II. THE MODEL 



Let us consider a laser pulse which first passes through a polarization beam splitter to separate two polarizations, 
e.g. horizontal and vertical, and then is sent through a fiber. Due to imperfections on the fiber, the two polarizations 
do not necessarily remain distinct, but they tend to become mixed. At the same time along the travel on the fiber the 
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FIG. 1: Polarization mixing and losses are modeled by three beam splitters characterized by real random variables 6 and <j> 
which have the probability distribution P(8,<f)). Input states enter the channel in modes 1 and 2 and the output is read in 
modes 7 and 8. 



pulse tends to be attenuated. This attenuation effect can be related to the polarization mixing effect. The situation 
can be modeled by three beam splitters as shown in figure {l}. 

In general the evolution operator of a beam splitter is given by the unitary transform 

[7(C) = e iabi - c * aib , (1) 

depending on a complex parameter £. Here a (a') and b (b') are annihilation (creation) operators for the two input 
radiation modes. It is easy to show that the Heisenberg evolution of the creation operators of the two modes of the 
input state is given by 



Ifi J u -e-iarsf sin|C| cos |C| J \ 

All the beam splitters in the channel are characterized by real parameters 9 and <j>. Labeling the creation operators of 
the input modes by a\ and a 2 and using equation ([2|), the creation operators of modes 3 and 4 in figure are given by 

ajj = cos 9a \ + sin 9a\ 

a\ = — sin 9a\ + cos 9a\ , (3) 

which shows that they are a mixture of modes 1 and 2, hence the polarization mixing effect. For example, if the mode 
1 (2) denotes the horizontal-polarized (vertical-polarized) photons, in mode 3 (4) photons with both horizontal and 
vertical polarization will appear. 

Then, each output of the first beam splitter goes through another beam splitter, characterized by a real parameter 
<j), and interacts with environment in the vacuum state. This part models the losses in the channel with a loss rate 
equal to sin 2 <f). 

The total action of these three beam splitters on the input state p\ t 2 is described by the map 

*fl,*(pi,a) = Tr e „„ [u 4fi (0)U 3 , 5 (<P)U h2 (9)( Ph2 ® |00) 3 , 4 <00|)[/t 2 (9)ul ^)ul , (4) 

where the environment modes are labeled by indices 5 and 6. To complete the definition of the channel we introduce 
correlation between polarization mixing (characterized by 9) and losses (characterized by <fr). This can be done by 
means of a joint probability distribution P(9, <fi) for the values that 9 and <j) can take. Thus, the output of the channel 
results 

$(p) = / dBd4>P{9,<t>)* g ,4>{p)- (5) 



We assume P(0, 4>) to be a Gaussian distribution of the following form: 
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where r — (6 — 9* , (f> — </>*)* and 7 is the correlation matrix: 

The parameter x € [0, 1] accounts for the degree of correlation between the two effects (polarization mixing and loss). 
For i = 0we have 

p{e,cj>)=p{6) P {4>) (8) 

where p{&) and p(<j>) are independent and identical probability distributions 

p(B) = -==e-^-, p(<t>) = -==e~^^ . (9) 

v 2,ira z v 2ira z 

They are characterized by mean values 9* and (j>* and variance a 2 . By increasing x, the distance between the 
probability distributions P{6, <fr) and p{9)p{4>) increases, 

\P(9, cf>) - p(6)p(<f>)\ = p{9)p{4>) J 1 - e -3**(«-«')(*-*> J . (10) 

This can be interpreted as increasing the correlation [13]. Since at x = 1 the maximum of the joint probability 
distribution is at 9 — 9* — <f> — <f>* we consider it as the most correlated case. 

The last remark which is necessary on the probability distribution © is that it is not a periodic function of 9 and 
4>, while the beam splitter operators are periodic in 9 € [— n, n] and <j) G [— 7r, 7t]. Hereafter we assume the variance a 
of the probability distribution to be sufficiently small (<§C 2ir) , so that boundary effects due to the periodicity become 
negligible. 

Our aim is to study the one-shot classical capacity for the channel ^ and analyse how it behaves with respect to 
the correlation parameter x. By definition, it is given by [HI]: 

C= sup xiHphPi)), (11) 
{pi,Pi} 

where the supremum is taken over all possible ensemble of input states, namelyjpi, pi} and \ is the Holevo information 

■(Pi)), (12) 



X (p) = s(*(y, Pip^j ^j-J2p* s ^ 1 



with S the von-Neumann entropy. 

In the following we are going to consider classical information over a continuous alphabet encoded in weak coherent 
states as well as classical information over a discrete alphabet encoded into qubit states. 



III. WEAK COHERENT INPUT 



Since coherent states turns out to be optimal for achieving the capacity in lossy bosonic channels Q, we consider 
them as input to the channel (O . A coherent state of complex amplitude a is expressed in the Fock basis by 



\a) = e—r^2-=\n). (13) 

We are here interested on weak (attenuated) coherent states which likely contain no photon or one photon, i.e. they 
are characterized by an amplitude |a| <C 1. To provide the two input modes of the channel in figure {l}, we consider 
horizontal and vertical-polarized photons so that the state of modes 1 and 2 in figure [T] can be written as 

\a) w JV(|00) 1>2 + a(|01 v )i, a + M> 1)2 )) (14) 

where N = , 1 represents the normalization factor. 

It is straightforward to show that the state (fH)) will be transformed to the following state after the action of the 
first beam splitter: 



JV[|00) +acos0(|Ol v ) + |1 H 0)) - asin0(|OlJ - |l v 0))] 3 ,4 
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cos%) 3 ,4 - A r [cos0|OO) + asm0(|Ol H ) - |l v 0))] 3 , 4 . (15) 



In the second step, the above state (|15|) interacts with modes 5 and 6 which are initially in the vacuum. Hence, we 
have 

iV[|0000) + a cos cos 0(| 01 v , 00) + |1 H 000» 

- acos0sin0(|OOOl v ) + |001 H 0)) 

- a sin cos 0(101^00) - |l v 000)) 

+ asin0sin0(|OOOl H )-|OOl v O))] 7 , 8>9 ,i O . (16) 

By introducing 

\f(a)) := |00) + a|01j-a|l v 0), (17) 
g{6,<t>) := 1 — (cos0 — sin 0) (cos — sin</>), (18) 

we can rewrite the state (fT6|) in the following form: 

\9) = N [ ff(0 ) 0)|OOOO)7,8,9,io + cos0cos0|a) 7 ,8|OO)9,io-cos0siii0|OO)7,8|Q;)9,io 

- sin0cos0|/(a)) 7i8 |OO) 9 ,io + sin0sin0|OO) 7 , 8 |/(a)) 9 ,io], (19) 

which shows that with some probability the input state \a) emerges at the output modes 7 and 8. However, there 
is a chance of loosing all the photons to the environment, and a chance that, though the photons do not go to the 
environment, the output state \f(a)) is different from the input one. 

To find the final state which represents the channel output we trace |^ r )(4 r | over modes 9 and 10 and then integrate 
over and as in ([5]), thus obtaining 

{ l + 2\a\ 2 (l-A-B) aD aD -o£ -o£ \ 

aD \a\ 2 A \a\ 2 A -\\ot\ 2 C ~\\a\ 2 C 

aD \a\ 2 A \a\ 2 A ~t\a\ 2 C -||a| 2 C 



$(\a)(a\) = N 2 



-a£ -\\a\ 2 C -\\a\ 2 C \a\ 2 B \a\ 2 B 

\ -a£ -\\a\ 2 C ~\\a\ 2 C \a\ 2 B \a\ 2 B J 



(20) 



where A, B,C,D and £ are real functions of 0*, 0*, a and x whose explicit form is given in Appendix. It is important 
to note that the order of \a\ is the same in the output state and the input state 

Once we have the output state of the channel we can then evaluate the Holevo information. We recall that the 
Holevo information for an ensemble & = {\ot}; p(a)} of input states corresponding to a continuous alphabet is given 
by: 

X (6) = S(p) - J d 2 ap(a)S (#(|a)<a|)) , (21) 
in which p is a probability distribution and p is the average output state 

p=$( / d 2 ap{a)\a)(a\ ) = / d 2 ap{a)<S> (\a) (a\) . (22) 



We restrict the ensemble of input states to that of weak coherent states (IT4")) with Gaussian distribution p(a): 

P{a)= 2^ e ~^- (23) 

The variance A 2 should be chosen in a way that the ensemble of input states only includes weak coherent states, i.e. 
A -c 1. We also know that and <fi are random variables with Gaussian probability distribution centered around 0* 
and cf>* with variance a. Therefore to have the channel close to an ideal transmitter, we assume narrow probability 
distributions (small a) centered around 0* = <f>* = 0. Under these assumptions C — £ w 0, therefore the output of 
the channel ([20|) is block diagonal with the following eigenvalues: 

JV 2 



A ll2 = — (1 + 2|«| 2 - 2\a\ 2 B ± ^(1 - 2\a\ 2 A) 2 + 8\a\ 2 V, 
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A 3 = 2N 2 \a\ 2 B, 
A 4 , 5 = 0. 

Using these eigenvalues the output entropy 



(24) 



Sm\a){a\))=-J2x i logX i 



(25) 



can be easily expressed analytically. Nevertheless, the integral in the second term of the Holevo information (j2~lj) 
should be computed numerically. 

For finding the first term of Holevo information we note that p(a) is even in a, therefore the average output state 
is block diagonal 



/1-2C(A + B) 



AC AC 
AC AC 



\ 



BC BC 
BC BC J 



(26) 



with the following eigenvalues: 



Ai = 1 + 2C(1- A-B), 

A 2 = 2 AC, 

A 3 = 2BC, 

A4.5 - 0, 



and 



C 



d 2 ap(a)- 



l + 2\a\ 2 ' 

Evaluating C numerically the first term of Holevo information in equation J21 



(27) 



(28) 



(29) 



can also be found. The final result for Holevo information (|21[) is plotted in figure ((2]) versus the correlation parameter 
x, for different values of a. It shows that the Holevo information for the ensemble of input states (f2"3"f is a decreasing 
function of correlation parameter x. For x — polarization mixing and loss happen independently, therefore there is 
the chance that just one of them happens. By increasing the correlation parameter the probability that both errors 
happen increases and therefore the rate of information transmission decreases. As it is clear from the figure ([2]), by 
decreasing a the Holevo quantity increases, this reflects the fact the narrow the probability distribution is, the less is 
the noise in the channel. As a goes to 0, A, T> — > 1 and B — > 0. In this case the output state gets closer to a pure 
state of the form of input state JIT 



IV. QUBIT INPUTS 



In this section we use qubit encoding in single photons. To this end we consider the two input modes 1 and 2 to 
only have one photon in total. Therefore the photon is either in mode 1 or in mode 2 and a general input pure state 
is given by: 

|V> = c |01> + Cl |10>, (30) 
with Co, Ci are complex coefficient such that |co| 2 + |ci| 2 = 1. By referring to the dual-rail coding 

|0>£ <-> |01} ll2 ^ |10)i, 2 . (31) 

it is clear that the state in ([30]) represent a general pure state for a single logical qubit. 
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FIG. 2: Holevo information versus correlation parameter x when information is encoded in weak coherent states. Different 
lines are for different values of a. 



Our aim is to find the capacity of channel ([5|) when we encode information in qubits. The upper bound of Holevo 
information (|12p is given by 

X<l-S(*(p*)), (32) 

in which p* is an input state with minimum output entropy [4[ . Using the concavity property of entropy, it has been 
shown that for finding the input state with minimum output entropy it is enough to search among the pure states 
[i[ . Therefore we consider the input ensemble that includes pure states of form ([30]) . Using equation (|2|) it is easy to 
show that 

17(0)|O1} = cos0|Ol) -sin0|lO), 

I7(0)| 10) = sin0|Ol) +cos0|lO), (33) 
therefore the action of the first beam splitter is a rotation by angle in the \0}l and plain, and the input state 

=co|0) L + ci|l) L , (34) 

will be transformed to 

\i>{6)) L = e- i6 °y\i,) L , (35) 

in which a y is the y Pauli operator. The dual-rail state at the output modes of the first beam splitter is then given by 

(co cos0 + Ci sin 0)101)34 — (cq sin0 — c\ 008 0)110)3,4. (36) 

This state interacts with modes 5 and 6 which are in the vacuum and will be transformed into 

(co cos0 + ci sin0)[sin</>|OOOl) + cos 0|O1OO)] 7 , 8 , 9,10 
- (c o sin0- Cicos0)[sin</>|OO1O) + cos0|1OOO)] 7 ^9,io- (37) 

If we take |2)i <-> 1 00) , then the output state of system and environment can be summarized as follows: 

sin0|2 s ) L |^(0) e ) L + cos0|^(0) s )L|2 e ) L , (38) 

where s and e are labels for system and environment respectively. Tracing over the environment and taking into 
account the correlation between the parameters of the beam splitters the final output state of the channel is: 

H\iP) L m = e\2) L {2\ + (1 - e)e- i6 ^y\i,) L {^ 6 '"y, (39) 
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FIG. 3: Classical capacity of the channel versus correlation parameter x when information is encoded into qubits. Different 
lines are for different values of a. 



where 

sin 2 (f)P{8, (j))d6d(j) = i ^1 -cos2</>*e - ^^ . (40) 

Since \2)l is orthogonal to the manifold of the input states, equation (|39p shows that the information is lost with 
probability e and the state is rotated by average angle 6* with probability (1 — e). Therefore the total map is a kind 
of a quantum erasure channel and its capacity is given by 

C=(l-e)C n , (41) 

where Cq is the capacity of the the channel 

n{p) = e- ie *' 7 y p e ie *' J y. (42) 
Since f2 is a rotation about y axis it is obvious that the following states remain invariant under its action: 

\^, 2 ) L = ^=(\0) L ±i\l) L ), (43) 

thus giving zero output entropy and minimizing the second term of Holevo information. 

It is easy to show that if these states are used with equal probability in the ensemble of input state then 

^(IV>iMV>il) + \m^)M) = \i. (44) 

Therefore this ensemble also maximizes the first term in Holevo information and Co = 1 can be achieved by using 
the above states to encode information. As a consequence the capacity of the channel $ results C = (1 — e). 

Figure ([3]) shows this capacity versus the correlation parameter x for different values of a when 9* =(/)*= 0. At 
x = the values of <fi which are close to are the most probable values, however by increasing x the higher values of 
<f> also have chance to happen. Therefore the probability of loosing information increases and as it is seen in figure ([3]) 
the capacity of the channel decreases. By decreasing a the probability distribution becomes more and more narrow 
and the chance for having values of <fi far from <p* = become less. Therefore the channel capacity increases by 
decreasing the variance of probability distribution which can be seen in figure f3|). 



8 



V. CONCLUSION 



We have provided a model of quantum channel with correlated errors of different kinds. While correlation in errors 
between different uses of channels have been studied before, in our model different types of correlated errors happen 
in each use of the channels on different degrees of freedom. This assumption can be justified by considering practical 
situations i. e. optical fibers where photon loss and polarization errors may be correlated. Actually our model 
describes photon losses depending on the degree of polarization mixedness. 

On the one hand, for classical information over a continuous alphabet we have considered weak coherent states and 
we have shown that the Holevo information decreases by increasing the correlation between polarization mixdness and 
losses. On the other hand, for classical information over a discrete alphabet we have considered single photons using 
dual rail representation. In this case we have shown that the channel is kind of quantum erasure channel where the 
probability of erasing information increases with the polarization mixdness. It is worth remarking that while often 
memory effects lead to an enhancement of the channel capacity, in the present model exactly the opposite happens. 
This is mainly due to the fact that the polarization mixing alone cannot be considered as a "true" error, i.e. it is not 
due to the action of an environment. 

In conclusion, we are confident that this work may pave the way for characterizing realistic quantum channels where 
different degrees of freedom and different effect are involved. 
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Appendix 

To evaluate the quantities A, B, C, T> and £ in ([2H|) we use the fact that the variance of the probability distribution 
P(9, 4>) is sufficiently smaller than 27T, hence we can extend the limits of the integrals to (—00, 00) treating them like 
standard Gaussian integrals, thus obtaining: 



A = J cos 2 6 cos 2 4>P(0, <j})d9d<j) 



= Ul + e ^ (cos 20* + cos 20*) + + e'^rl ) cos 20* cos 20* j , 

B = J sin 2 0cos 2 0P(0,0)d0d0 

= - fl - e'^ 1 (cos 29* - cos 20*) - + e'^rl ) cos 20* cos 20*^j , 

C = J sin20cos 2 0P(0, 4>)d9d(j) = ~ ^2e~^ - (e~^ + e^^rl) cos0*^j sin20*, 

V = fcos0cos0P(0,0)d0# = i (e~& + e"^) cos 0* cos 0*, 

£ = /sin0cos0P(0,0)d0d0 = - (e~^ +e~TT^ sin0*cos0*. (45) 
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